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CN ' Spin waves in (0, vr) and (0, vr, vr) ordered spin-density-wave (SDW) states of the 

, ^ . t — t' Hubbard model are investigated at finite doping. In tlie presence of small 

^O I t', these composite ferro-antiferromagnetic (F-AF) states are found to be strongly 

stabilized at finite hole doping due to enhanced carrier-induced ferromagnetic spin 

^ ■ couplings as in metallic ferromagnets. Anisotropic spin-wave velocities, spin-wave 



energy scale of around 200meV, reduced magnetic moment, and rapid suppression of 
magnetic order with electron doping x (corresponding to F substitution of O atoms 



'^ ■ in LaOi_xFxFeAs or Ni substitution of Fe atoms in BaFe2_xNixAs2) obtained in this 

O . model are in agreement with observed magnetic properties of doped iron pnictides. 



PACS numbers: 75.30.Ds,71.27.+a,75.10.Lp,71.10.Fd 



I. INTRODUCTION 

Following the recent discovery of superconductivity^"- in doped iron pnictides such as 
ROi_xFxFeAs (R = La, Ce, Nd, Sm, Gd), there has been a strong interest in their magnetic 
properties as well. Single-crystal neutron scattering studies of their parent compounds have 
indicated a commensurate magnetic ordering with iron moments ordered ferromagnetically 
in the b direction and antiferromagnetically in the a and c directions.- Inelastic neutron 
scattering measurements in AFe2As2 (A = Ca, Ba, Sr) yield sharp spin-wave excitations on 
an energy scale ~ 200meV.-"- 

AU known compounds in these classes contain FeAs layers and exhibit a common phase 
diagram, with parent compounds exhibiting a magnetically ordered SDW state below 
T/v ~ 200 K, and onset of superconductivity following the suppression of long-range mag- 
netic order by electron doping or pressure. Contrast with cuprates has naturally followed 
in order to gain insight from the significant progress made in understanding superconduc- 
tivity and magnetism in correlated electron systems.- While onset of superconductivity at 
finite doping is a common feature, in contrast to the Mott insulating nature of cuprates, the 
pnictides appear to be commensurate SDW systems in the intermediate coupling regime. 
Appreciable hybridisation between Fe 3d orbitals and As 3p orbitals possibly reduces the 
effective correlation term U as compared to cuprates.-™ Comparison with X-ray photoemis- 
sion spectra of the calculated density of states (DOS) within the LDA + DMFT approach 
supports the physical picture of a multi-band metal with intermediate correlations.—*^ 

In cuprates, the intense interest in the nature of magnetic excitations in the quantum 
antiferromagnet, their coupling with hole motion and scrambling of local AF order, strongly 
incoherent hole spectral function, and possibility of pairing interaction mediated by ex- 
change of local magnetic excitations have contributed significantly to obtaining insight and 
understanding of their magnetic and electronic properties. 

Therefore, an investigation of magnetic excitations in the (0, vr) and (0, vr, vr) ordered 
SDW states within a minimal itinerant electron model should be of interest, particularly 
within an approach which is valid in the full range of interaction strength including the 
relevant weak and intermediate coupling regimes. In this paper, we will therefore consider 
the t — t' Hubbard model, and obtain the magnon energies in the F-AF states in two and 
three dimensions, focussing especially on the role of finite t' and doping in stabilizing the 



F-AF order. 

Spin waves in the (0, vr) and (0, n, n) ordered states were investigated earlier in the un- 
doped t — t' Hubbard model,— and have been recently investigated within Ji — J2 and 
Jia — Jib — -h Heisenberg models on a square lattice.— ^^ Reduced magnetic moments and 
suppression of magnetic ordering have been investigated in terms of associated magnetic 
frustration effect in such models.— ^^^ Spin excitations in SDW states have been investigated 
within two-band models involving excitonic instability of nested electron-like and hole-like 
Fermi pockets.— Doping dependence of spin fluctuations and electron correlations have been 
theoretically investigated within an effective five band Hubbard model for iron pnictides us- 
ing the FLEX approximation.— 

II. THE t - t' HUBBARD MODEL 

We consider the t — t' Hubbard model on square and simple cubic lattices, with hopping 
terms t and t' between nearest-neighbour (NN) and next-nearest-neighbour (NNN) pairs of 
sites, respectively: 

NN NNN 

H = -t^ al^ai+s,a -t'^ a\a"'i+K,a + f^ XI ^^^^4 • (1) 

The observed asymmetry in the antiferromagnetism of cuprates with respect to hole and 
electron doping,— 1^ existence of metallic antiferromagnetism at half filling,— 1^ and corre- 
lated motion of electrons in metallic ferromagnets in fee and bcc lattices,—"— exemplify few 
physical systems which have been investigated in terms of the above model. 

The two-sublattice representation of the AF state (corresponding to ordering wave- vector 
(vr, tt) in two dimensions) conveniently allows for investigation of spin waves, quantum cor- 
rections, Neel temperature, hole dynamics etc. In analogy, we will consider F-AF SDW 
states with Q = (0, vr) and (0, vr, vr) involving ferromagnetic spin ordering in one direction 
and antiferromagnetic spin ordering in the remaining direction(s). The self-consistent-field 
(Hartree-Fock) approximation provides a convenient basis in which many-body correlations 
effects can then be systematically incorporated. The two-sublattice structure for the F-AF 
SDW states in two and three dimensions are given below. 



A. Q = (0, vr) ordered state 



In this state, spins are ordered ferromagnetically along the x direction and antiferromag- 
netically along the y direction. The NN hopping terms in the x direction connect sites of 
same sublattice, while NN hopping terms in the y direction and NNN hopping terms con- 
nect sites of opposite sublattices. The HF Hamiltonian matrix in the two-sublattice basis 
therefore takes the form: 



^SF(k) 



-aA 



el + el 



el + e; aA + el 



^kl 



-a A r]k 
^k o-A 



for electron spin a, where ^k = ^k ~ ~2t cos k^ and r^k 



e^ 4- e' 
e^r \ fct 



4t' COS kx COS ky, and the self-consistent exchange field is given by 2A 
the sublattice magnetization m. 

The SDW state quasiparticle band energies: 



(2) 



= — 2t cos ky — 
mil in terms of 



Eit^ =vi± ./ A2 + r^l 



^ko- 



(3) 



corresponding to the two spins (c =t, J,) and the two (±) bands, and the quasiparticle 
amplitudes Cko- and foko- on the two sublattices A and B are given by: 

A 



'^kte ~ ^k^e 



''HQ 



M® 



'■kfa 



"k4.g 



"k^e ~ "kte 




(4) 



These relationships follow from the spin-sublattice and particle-hole symmetry in the AF 
state. The above two expressions provide the majority and minority fermionic densities. On 
the A-sublattice, the majority density is of spin f (J,) states in the lower (upper) band. 



B. Q = (0, TT, vr) ordered state 



In this state, spins are ferromagnetically ordered along the x directions and antiferro- 
magnetically ordered along the y and z directions. The NN hopping terms in the y-z plane 
connect sites of opposite sublattices, while those in the x direction connect sites of the same 
sublattice. Similarly, NNN hopping terms in the y-z plane connect sites of the same sub- 



lattice while those in the z-x and x-y planes connect sites of opposite sublattices. The HF 
Hamiltonian matrix therefore takes the form: 



^SF(k) 



where 



-a A + el + e'J' 



et + el + e' ^" + e 






/ xy 
k 



V I z \ I zx , I xy 

Ck + ^k + Ck + ek 



aA + e? + e'/^ 



^kl + 



-crA r/k 
r/k o-A 



'7k — ^k + ^k^^ = ~"2t COS A;2; — At' cos /cj, cos kz 

'7k = ^k + ^k + ^k^^ + ^''k^ = [~2t — At' COS A;^] (cos /cj; + cos kz) 

The quasiparticle band energies and amplitudes are again as in Eqs. (3,4). 

III. MAGNON PROPAGATOR 



(5) 



(6) 



Magnon excitations in the F-AF state, corresponding to transverse spin fluctuations about 
the ordering direction (assumed z), are obtained from the time-ordered propagator of the 
transverse spin operators 5*^" and S^ at sites i and j: 

X-+(q,a;) = /" rft J^e*'^^*-*^-^'!-^^-'-^) {^>o\T[S-{t)Sl{t')]\mG) • (7) 

In the random phase approximation (RPA), the magnon propagator: 

-+ ( \ [x°(q,^)] (^s 

XHPA(q.-)=i_j,|^„(q,,)| (8) 

in terms of the bare particle-hole propagator [x°(q, w)], evaluated by integrating out the 
fermions in the spontaneously-broken-symmetry F-AF state. In the insulating state, involv- 
ing only interband particle-hole processes, the bare propagator is given by: 

Okte ^kte ak-q4.® ^k-ql® 

akte^ktectk-q4.ffi&k-q4e ^kte^k-q^e 

'^kte'^k-qie Okte^kteOk-q^e^k-qie 

akte^ktffiCtk-q^e^k-q^e ^kt®^k-q^e 

in terms of the quasiparticle amplitudes and energies. In the antiferromagnetic metallic 
(AFM) state, additional intraband processes involving particle-hole excitations from the 
same band also contribute.—"^ Evaluation of [x°(q, w)] and the RPA-level magnon propa- 
gator x~+(q, cj) in the strong coupling limit is described in the next section. 



[X°(q,c.)] = Y. 

k 



'^kte'^k-q^.g 



1 



Et 



^k-q - ^ 



17] 



17] 



(9) 



IV. INSULATING (0, tt, vr) STATE IN THE STRONG COUPLING LIMIT 



In this section we consider the analytically simple strong coupling limit, and evaluate the 
magnon propagator in the insulating (0, vr, vr) state at half filling, illustrating the competi- 
tion between NN (J = At^ /U) and NNN (J' = 4t'^/U) AF spin couplings and consequent 
instability of this state as J' drops below J/4. 

In the strong coupling limit, the majority and minority quasiparticle densities in Eq. 
(4) reduce to approximately 1 — ri^/AA"^ and r7^/4A^, respectively. Similarly expanding the 
energy denominators involving the quasiparticle band energies in powers of t/A, t'/A, and 
uj/A, and systematically retaining terms up to order t^/A^ and t ^/A^ in the bare particle- 
hole propagator, we obtain for the RPA level magnon propagator:— 



[x ^(q,^)] 



■-(-) 



A — — 

i3n 



B^ 



vAq 



q ' 2j 

where the magnon propagator matrix elements: 
2J'\ I { ,_ . 2J' 



1 



1 



w — Wq + ir] uj + UJc^ — ir] 



(10) 



A, 



Sn 



J 

2J' 



1 — COS Qx 



J 



1 — cosgy cosg^) 



cos Qa 



(cos qy + cos Qz 



J 'J 2 
and the magnon-mode energies are obtained as: 

2 



(11) 



2j) 



Ai-Bi 



2J'\ \ { , . 2J' , 

1 + - -2 (l-cosg.) + -(l 



cos qy cos qz 



2J' 



cos qx 



cos qy + cos qz 



J ' 'V \ 2 

In the long wavelength limit, the magnon energy reduces to 



m 



2 1 / 2J' 

1 + ^1 [aqi + q'y + q 



.^], 



(12) 



(13) 



2 V J , 

where the coefficient a = (^j — l) of the g^ term becomes negative for J' /J < 1/4, signaling 
the instability of the F-AF phase at J' /J = 1/4. Anisotropic spin wave velocities naturally 
follow from the different q^ coefficients in the F and AF directions. 

The above instability can also be seen from energy considerations. The classical energy 
per spin for the two orderings are: Ea_p = —QJ + 12 J' and Ep-^F = —2 J — 4 J', so that the 
F-AF state becomes energetically favourable for J' > J/4. In three dimensions, the colinear 



1 1 1 1 1 1 1 1 r 




F-AF (0,7r,7r)- 
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FIG. 1: The spin-fluctuation correction to sublattice magnetization in the AF and the F-AF phases. 

Q = (0, TT, 7r) state is stable even at the classical level, unlike the degeneracy present in the 
d = 2 case at this level. 

As an illustration of quantum corrections beyond the HF approximation, the spin- 
fluctuation correction ^tt^sf to the sublattice magnetization in the F-AF phase, which can 
be evaluated from the magnon propagator in terms of the transverse spin correlations,— is 
shown in Fig. 1. The correction in the AF phase is also shown for comparison. Near the 
transition point J'/J = 1/4, the correction in the F-AF phase is seen to be nearly half 
of that in the AF phase, indicating greater robustness of the F-AF phase with respect to 
quantum spin fluctuations. The spin-fluctuation correction in both phases approaches 1 (the 
HF value of sublattice magnetization) only very close to the critical value J'/J = 1/4. This 
implies that (up to first order) m vanishes only very close to J'/J = 1/4, so that the extent 
of the spin-disordered phase is quite narrow. This is unlike the d = 2 case, where the AF 
order is lost at J'/J ~ 0.37, well before the F-AF state appears at J'/J ^ 0.5. 

V. STABILIZATION OF THE HOLE-DOPED F-AF STATE 



As shown by the strong-coupling analysis of the previous section, AF NNN spin couplings 
ji _ 4^'2y'f/ generated by t' stabilize the undoped F-AF state (for t' /t > l/-\/2 in two and 
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t' jt > 1/2 in three dimensions). So how is the F-AF state affected by other physical processes 
which might also generate effective NNN spin couplings? 

AF NNN spin couplings are generated even in the absence of t', as in the t — U model at 
finite doping. In fact, it is these effectively frustrating spin couplings, generated by intra- 
band particle-hole processes, which are responsible for destabilizing the AF state at any 
finite doping, as observed in hole doped cuprates.-2i However, as the F-AF state is actually 
stabilized rather than being frustrated by the AF NNN spin couplings, these doping-induced 
spin couplings should actually favour the F-AF state by supplementing the t'-induced spin 
couplings. 

More importantly, carrier-induced F NN spin couplings responsible for metallic ferro- 
magnetism become increasingly important at finite doping. This is especially so in presence 
of small t', which can cause strongly peaked electronic spectral distribution due to band 
dispersion saddle points (Vck = 0), which strongly enhance band ferromagnetism by in- 
creasing the delocalization contribution ~ (V ek) to spin stiffness while strongly suppress- 
ing the correlation-induced exchange contribution ~ ((Vek)^/f/) due to correlated motion 
of electrons. ^^'^^ The F-AF state at finite doping is therefore expected to be stabilized at 
even lower t' values. In this section we will show that indeed finite hole doping and small t' 
strongly stabilizes the F-AF state. 

Fig. 2 show the spin wave energy along symmetry directions in the Brillouin zone for the 
doped F-AF state, with orderings (0, vr) and (0, vr, vr) as considered earlier. The variation of 
wave-vector q follows the sequence (0, 0) — t- (0, vr) — )■ (tt, tt) — ;■ (0, 0) -^ (vr, 0) — )■ (vr, vr) in 2d 
and (0, 0, 0) -^ (0, vr, tt) -^ (0, tt, 0) -^ (vr, vr, 0) -^ (vr, 0, 0) -^ (0, 0, 0) in U. In order to focus 
on doping dependence, the SDW bands were kept unchanged with fixed A/t = 4. The rapid 
crossover from negative to positive energy magnon modes shows a strong stabilization of the 
F-AF state upon hole doping. This stabilization occurs for much smaller t' values compared 
to the critical values required in the undoped F-AF state. A finite t' is quite realistic in view 
of the hybridization between Fe and As orbitals. 

The anisotropic spin wave velocities, evident from the different slopes in the F and AF 
directions in Fig. 2(a), can be understood readily in terms of the independent origin (de- 
localization and exchange) of the effective F and AF spin couplings, resulting in different 
coefficients of g^ and g^, as in Eq. (13). Furthermore, for t ~ 200meV and U in the in- 
termediate coupling regime, the calculated spin wave energy scale of around 200meV is as 
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FIG. 2: Spin wave energy along symmetry directions in the Brillouin zone for the F-AF states 
with (a) (vr, 0) and (b) (0, vr, vr) orderings. The rapid crossover from negative to positive energy 
modes, especially in the F direction F-M', shows a strong stabilization of the F-AF states upon 
hole doping y. Here t' /t = 0.5 in (a) and 0.3 in (b). 

observed in neutron scattering measurements of iron pnictides. 

Are the doping-induced spin couplings sufficient to stabilize the F-AF state without 
any finite t'7 We find that in the absence of t' the F-AF state is not stabilized for any 
doping. While AF and F orderings do get separately stabilized at low and high hole dopings, 
respectively, as indicated by the spin-wave dispersion along AF and F directions, both are 
not simultaneously stabilized. Thus, small t' plays a crucial role in stabilizing the doped F- 
AF state with respect to transverse spin fiuctuations in the entire Brillouin zone. A possible 
link between orthorhombic distortion and an effective NNN hopping t' would then explain 
why this distortion appears necessary for magnetic ordering to be stabilized in iron pnictides. 

As the maximal sensitivity of magnon modes to doping is found along the ferromagnetic 
direction (F-M'), we have examined their doping dependence at an intermediate wavevector 
q = (vr/2, 0, 0). A pre-doped level with hole doping yo was taken to represent the partially- 
filled band of the parent compound, with electron doping x = yo — y defined as the reduction 
in hole doping from this level (representing F~ substitution of 0^~ in the doped pnictides). 
Fig. 3 shows the behaviour of the magnon mode energy with electron doping. The rapid 
magnon energy suppression and crossover to negative energy modes indicating destabiliza- 
tion of the F-AF state provides an understanding of the observed rapid suppression of 
magnetic ordering temperature in iron pnictides with electron doping (due to F substitution 



10 





120 




100 


> 


80 
60 


s" 


40 




20 









FIG. 3: Electron doping (x) dependence of the magnon energy at wavevector q = (7r/2,0, 0), 
showing rapid suppression of magnon mode energy along the ferromagnetic ordering direction and 
destabilization of the F-AF state. 

of O atoms in LaOi_xFxFeAs or Ni substitution of Fe atoms in BaFe2_xNixAs2).^° Here, the 
interaction strength was fixed at U = 8t in the intermediate couphng regime, and yo = 0.35. 
The observed reduced magnetic moment in pnictides can be understood in terms of 
depletion of the predominantly magnetic states. Characteristic of the SDW state, electronic 
states at the top of the band are significantly more magnetic than states deeper in the 
band, especially in the weak-coupling limit. Therefore hole doping of these states rapidly 
diminishes the local magnetic moment m = 2A/U, which is reduced to ~ 0.4 at x = in 
Fig. 3. While local moments will get enhanced on filling up these empty magnetic states 
by electron doping (F substitution of O), the rapidly diminished carrier-induced F NN spin 
couplings are then unable to sustain the F-AF ordering. This highlights the two distinctly 
different mechanisms behind the observed reduced magnetic moment and the suppression of 
magnetic ordering in iron pnictides. Enhanced spin-fiuctuation quantum correction in the 
vicinity of the magnetic instability point (as in Fig. 1) is also possibly important in reducing 
the magnetic moment. 

VI. CONCLUSIONS 



Spin waves in the (0, vr) and (0, vr, vr) ordered F-AF SDW states of the t — f Hubbard 
model were investigated at finite doping in the intermediate coupling regime. Including both 
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inter-band and intra-band processes in the particle-hole propagator to obtain the spin-wave 
propagator which incorporates effective spin couplings, spin-wave dispersion was obtained 
along different symmetry directions in the Brillouin zone. The F-AF state was found to be 
strongly stabilized at finite hole doping, as evidenced by the rapid crossover from negative 
to positive spin-wave energies, most noticeably in the F direction. 

This stabilization was ascribed mainly to carrier-induced ferromagnetic spin couplings 
(along the F chains) as in metallic ferromagnets, and was found to be strongly enhanced 
in presence of finite t' due to band structure saddle point effect. This calculated doping 
behaviour is in agreement with the observed rapid suppression of the magnetic order in 
doped pnictides such as LaOi_xFxFeAs on electron doping arising from F substitution of 
O atoms. The doping and t' dependence of spin waves obtained within the minimal t — t' 
Hubbard model accounts for many of the observed magnetic properties of iron pnictides, 
including anisotropic spin wave velocities, spin wave energy scale, reduced magnetic moment, 
and rapid suppresion of magnetic ordering with electron doping. The spin-wave velocity 
anisotropy is strongly enhanced on electron doping near the instability point. 

Quantum corrections to spin waves beyond RPA should be of interest in order to in- 
corporate the correlated motion of electrons, as investigated in the metallic ferromagnetic 
state recently,— i^^ and the insulating AF state earlier,™ incorporating correlation-induced 
self-energy and vertex corrections within a systematic inverse-degeneracy (l/A/") expansion 
scheme which preserves spin-rotation symmetry and hence the Goldstone mode at each 
order. 

In an orbitally degenerate metallic ferromagnet, the inter-orbital Coulomb interaction 
(Hund's coupling) J was recently shown to strongly suppress quantum corrections,^^ partic- 
ularly for large orbital degeneracy A/", with the quantum correction magnitude determined 
by an effective quantum parameter J J^j^^J^, j / J21 • The renormalized spin stiffness with re- 
alistic parameters was obtained in close agreement with the measured spin stiffness for bcc 
Fe. In iron pnictides, the reduced effective A/" due to hybridization between (magnetic) Fe 
and (non-magnetic) As orbitals should enhance quantum corrections and hence suppress 
(F-direction) spin wave energies and the ordering temperature. These considerations are 
relevant in view of recent magnetic form factor studies indicating that multiple d orbitals of 
iron atoms are occupied.— 

The magnon energy scale also determines finite temperature magnetic properties. 
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Magnon thermal excitation yields the fall off of magnetization with temperature and, within 
a self-consistently renormalized spin fluctuation theory, the magnetic ordering temperature. 
Thus, the rapid suppression of magnon energy with electron doping obtained in this paper 
accounts for the observed reduction of magnetic ordering temperature on F substitution of 
O atoms in LaOi_xFxFeAs. Furthermore, spatially anisotropic magnetic couplings, with 
small ratio {J_l/J\\ = r ^ 1) of interlayer to planar magnetic couplings, reduces the order- 
ing temperature to ~ T^'^°/ln(l/r) compared to the isotropic case, which possibly accounts 
for the low ordering temperature observed in doped pnictides (~ 200K), as in the layered 
cuprate antiferromagnet La2Cu04 where J ~ 1500 K and magnon energy scale ~ 250meV, 
but the Neel temperature is only about 400K. 
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